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ELECTROSTATIC ION CYCLOTRON INSTABILITY DUE TO A NONUNIFORM 
ELECTRIC FIELD PERPENDICULAR TO THE EXTERNAL MAGNETIC FIELD 


I. Introduction 

The ion cyclotron instability has been of importance to both space and 
laboratory plasmas. In most previous studies, field aligned currents or 
ion beams have been cited as the driving mechanism*. Recently some 
laboratory experiments^ have reported ion cyclotron instability in 
circumstances in which neither of the above mechanisms provides a 
satisfactory free energy source, and the existence of a two dimensional 
electric field is conjectured to play a role. Further, unstable ion 
cyclotron waves have been reported in connection with double layers'* and 
electrostatic shocks 1 *. In all these cases a localized electric field 
perpendicular to the external magnetic field is an intrinsic feature of the 
equilibrium. The purpose of this letter is to report a new kinetic ion 
cyclotron instability driven by a nonuniform electric field perpendicular 
to the external magnetic field which could be of importance to the 
situations referred above^”^. 


II. Theory and Results 

Using physical arguments we first establish that an equilibrium 

characterized by an external uniform magnetic field B Q (in the z direction) 

perpendicular to an external nonuniform electric field E (x) (in the x 

o 

direction) is unstable to the electrostatic ion cyclotron waves. The local 
dispersion relation of the electrostatic ion cyclotron waves^ is given by, 


1 + 1 + 1 + 1 !-nrTr' z 'iir*r' - ° <» 


2 

where (kX.) « 1 is assumed. Also t 3 T./T , v 4 Q are the ion and 

i i e i> e 

electron thermal velocities, = eB Q /mjC is the ion gyro 
radius, kj and k^ are the parallel and perpendicular components of the wave 

vector (with respect to B Q ), T (b) 3 I (b)exp(-b), 

2 2 n n — — 

b_ » k^ p^ /2, p^ 3 and *n are the modified Bessel's functions. 

Setting k ♦ 0 one obtains the dispersion relation for the Bernstein modes 

(BM), 


D RM (w,k) 3 I 


( 2 ) 
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and if (Wv g ) « k H « (w - nP^/v^ then one gets the corresponding 
expression for the ion cyclotron modes (IC), 


:a 2 r 


DjcU.k) 


1 + T ' 


0 . 


(3) 

local 


L 2 2 2 ' 
n>0 m - n P 

For the time being we neglect all the natural dampings. In the 
approximation the energy density of the electrostatic waves is proportional 
to 3(u>D)/3ui «* ii3D/3aj. Thus the energy density U, of both the BM & IC waves 
is given by 


U « ui(3D/3to) 


4r n Q , 

{r U) —- u>"a(o>), a > 0. 
n>0 (a> - n fl“)“ 


(4) 


Clearly U is positive definite and thus the waves are positive energy 
waves. 

We now introduce a uniform electric field E Q in the x direction. This 
initiates a drift of magnitude Vg = cE/B in the v direction. Thus, other 
than a Doppler shift (i.e., u> ♦ * to - k^Vg) there is no change in the 

dispersion relations given in (2) and (3). The energy density of the 
Doppler shifted waves V' , is given by, 

U" * a) o)jo(o)j) . (5) 

While U is positive definite IT can be negative 

provided 0 < co < k y v p thereby giving the waves a negative energy density. 

Clearly a uniform electric field perpendicular to the external 
magnetic field can convert the positive energy waves into negative energy 
waves. However, as long as the electric field is uniform, a transformation 
to another frame moving with a velocity Vp will enable the waves to get 
back the positive energy character. This is no longer possible once the 
electric field is nonuniform. For example consider a specific model for 
the external electric field, 

E o > - i < * ‘ 7 

E q (x) = j 

0, otherwise 


2 








As the main objective of this letter is to demonstrate and discuss the new 
instability we have chosen a somewhat idealized electric field profile for 


relative ease in computation and transparency of the physics involved. 
Since Z is localized in x over a distance L (L > , one finds that a 

negative wave energy region I, is surrounded by positive wave energy 
regions II (See Figure 1). A nonlocal wave packet can couple these regions 
and a flow of energy from region I to region II can enable the mode to 
grow. This gives rise to the instability. This is essentially a nonlocal 
instability with no local limit. The negative energy wave framework has 
been previously used^ to explain other plasma instabilities at a local 
level. Ours is a nonlocal instability similar in concept with Lau^ at. al 
who apply this concept to an astrophysical problem. 

Since the x direction is nonuniform we use nonlocal differential 
equations 5 in the regions I and II to analyze this case. Dissipation 
processes such as Landau damping are now included. 

a 2 2 

REGION - 1: + k*\ $,(£) - 0, (6) 

ar 

a 2 2 

REGION - II: j — 2 + T }* TT (E) = 0, (7) 

H 

2 

where 5 = * Qj/A^ and 


b>. to. - nfl 

,___Wf —_; 

n ' |k ( |v t > ' |k ( | Vi 




Ik,|v ' 1 |k, |v 


A I = ” 2 n F n ' |k # Iv^ |k B |v ± ^ 


r 3 _ . ,2 2 

r = —r , b * k p. i . 

n 3b n’ y M i 


( 8 ) 


(9) 


is identical to if is replaced by u . The solutions in region I 
is, 
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$l(5) = fl^cos^E) 


(10) 


while in region -II it is, 

$^(0 = $ oI][ exp(iic II C + 16) » (11) 

where and 6 are constants. Demanding that the logarithmic 

derivative of the solutions be continuous at E = L/2p^ we obtain the 
nonlocal dispersion relation, 

-Kj tan(i<r I /2e) * iie^, (12) 

where e = p^/L. Given the parameters, we use a numerical root finder to 
solve (12) for the complex frequency <d. As a boundary condition we demand 
that the solution vanish at », i.e., InKic^) > 0. It should be remarked 
that since the electric field has a sharp profile our treatment at the 
second order differential equation level can be improved by employing an 
integral equation formulation. However, such a formulation is 

mathematically cumbersome and does not provide a simple and clear 
interpretation of the physics involved. Thus, for the time being, we 
stress simplicity and clarity over greater accuracy which will be provided 
in the future. 

Figure 2 is a plot of the real and imaginary parts of the 
frequency, u> r /Q and y/Q, against b for t = 0. , V £ = 2. 9v^, e ■ 0.3, 
u = k u/ky * 0.001 and y * m j/ m e = 1837. The instability exists only for a 
very narrow band of k y thereby making it very coherent. The real part of 
the frequency is large (<u - 1.60fl^) compared to the first harmonic. Also 
the instability peaks for ko^ *• 0(1). These features are In keeping with 
the observations of Mozer 4 et. al. 

Figure 3 is a plot of y/0,^ against u for t * 0, 0.5 and 1.0 and b * 
0.57 while other parameters are similar to the Figure 2. Clearly this 
instability prefers small u (perpendicular propagation) .again in keeping 
with Mozer et. al\ For t * 0, u can get as high as 0.15 but for 
higher x the value of u is even smaller. For u - 0.01 larger x implies 

4 



more electron Landau damping which reduces the growth rate. 


Figure 4 is a plot of go against e for u = 10 - ^ and other parameters 
similar to Figure 2. Our theory at the second order differential equation 
level is valid for e < 1 (i.e., p i < L) . For p^. > L situations an 
integral equation formulation will become necessary. For most of the 
laboratory experiments L *• (2 to SJpj , for which our theory holds. Figure 
3 indicates that the instability is more severe for smaller L (i.e. , for 
strongly localized electric fields) and disappears for L > 70^, for the 
given parameter range. 

Figure 5 is a plot of against V £ (* Vg/v^ for b = 0.57, 

t = 1.0, e = 0.3 and u = 10”^. The instability exists for V - between 2.3 

_ E , 

and 3.2 peaking at around 3. The peak value of V p in both Mozer et al. 
and Merlino et al. ^ is well above the minimum necessary for this 
instability. Further, is necessary for the onset. 

III. Conclusions 

In this letter we give a new mechanism that can explain the occurance 
of the short wavelength (kp^ - 1) turbulence around the ion cyclotron 
frequency in the presence of a nonuniform electric field perpendicular to 
the external magnetic field. This instability is driven by the coupling of 
the positive and negative energy ion cyclotron waves. Ion cyclotron 
turbulence has been observed associated with shocks^ in the magnetosphere 
and with double layers^ both in laboratory and space plasmas, where a 
strongly localized electric field perpendicular to the magnetic field is 
present. Thus the instability mechanism described in this letter could 
play a significant role in such situations. We would also like to note 
that the mechanism of the instability given here is general and can be 
applied to other instabilities as well. For example, a study of lower 
hybrid wave excitation by this mechanism is In progress and will be the 
subject of a future article 
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A plot of against V^. Other parameters are similar to 

the previous figures. 
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